Slow flows of yield stress fluids: 
complex spatio-temporal behaviour within a simple elasto-plastic model 
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A minimal athermal model for the flow of dense disordered materials is proposed, based on two 
generic ingredients: local plastic events occuring above a microscopic yield stress, and the non-local 
elastic release of the stress these events induce in the material. A complex spatio-temporal rheo- 
logical behaviour results, with features in line with recent experimental observations. At low shear 
rates, macroscopic flow actually originates from collective correlated bursts of plastic events, taking 
place in dynamically generated fragile zones. The related correlation length diverges algebraically 
at small shear rates. In confined geometries bursts occur preferentially close to the walls yielding 
an intermittent form of flow localization. 

PACS numbers: 83.50.-v, 81.40.Lm,62.20.Fe 
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Disordered dense systems often exhibit a peculiar flow- 
ing behaviour which strongly departs from the academic 
Newtonian description. A non linear rheology, with a 
shear rate dependence of the viscosity, or the existence of 
a yield stress are typical features of such systems. More- 
over, it has been recognized recently that these global 
characteristics are, in most cases, associated with a pe- 
culiar spatial behaviour, in the form of heterogeneous 
flow behaviour, where a frozen region coexists with a 
flowing one (the so-called "shear band"). A striking re- 
mark is that such generic behaviours are observed in a 
wide class of experimental systems, with very different 
length/time/interaction scales, such as foams gran- 
ular systems 0, , emulsions H, H, 1^ i colloidal glasses 
0, polymers, but also in simulations of granular systems, 
foams and model glasses lOj 11]. These generic fea- 
tures suggest an underlying common scenario for the flow 
properties, and has motivated various macroscopic phe- 
nomenological approaches (see eg refs. cited in (3] and 
0). However a consistent framework linking the global 
rheology to the local microscopic dynamics is still lacking, 
although some prog ress in this direction has been made 
in recent years |l2L lisl . In particular studies have 
put forward the role of local plastic rearrangements in 
the global flow behaviour |l11,.15, 16] . Such an idea actu- 
ally goes back to the Princen model for the deformation 
of foams : flow occurs via a succession of reversible 
elastic deformations and irreversible plastic events ("Tl" 
events in foams), associated with the existence of a lo- 
cal yield stress. However, if the corresponding physical 
picture seems a priori quite clear, a gap still persists be- 
tween this simple microscopic scenario and the complex 
spatio-temporal organization responsible for the rheology 
of these materials at finite shear rates. 



In this letter, we propose a simple mesoscopic model, 
constructed on the basis of two minimal and generic in- 
gredients : localized plastic events associated with a mi- 
croscopic yield stress, and the resulting elastic relaxation 
of the stress over the system. We then show that the 
simplicity of the description contrasts with the complex 
rheological behaviour deriving from it. In particular we 
find the global rheology to be associated with a complex 
spatio-temporal organization which builds up as the sys- 
tem is sheared steadily, with an intermittent behaviour 
corresponding to "bursts" of correlated events, the typ- 
ical size of which diverges at small shear rate. We ar- 
gue that in its present simple form our model seems to 
capture many observed experimental features and thus 
stands as a promising starting point for the elaboration 
of a generic scenario for the slow flow of yield stress fluids. 

Let us now precise the ingredients of our approach, 
which we implement here in the simplifled frame of a 
2D scalar approach, focusing only on the simple shear 
components of the stress and strain. We consider a two 
dimensional material to which a global shear rate 7 is 
applied (corresponding to a z dependent displacement in 
the X direction). The material is described at a coarse 
grained level, intermediate between the microscopic (par- 
ticle) and macroscopic scale. The quantity of interest is 
the xz component of the time dependent local shear stress 
a{x, z; t). First, without entering into details at this level, 
a few basic rules are stated : (i) below a (locally defined) 
yield stress cry, the system responds elastically to the 
imposed deformation; (ii) above ay, plastic events may 
occur in the system (along laws discussed in the follow- 
ing); (iii) plastic events take the form of a localized shear 
strain; (iv) such a plastic event induces a long range elas- 
tic perturbation of the shear stress field in the material. 
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A few remarks can be done at this level. First, although 
the notion of individual events is quite intuitive, in partic- 
ular in foams, it has been evidenced unambiguously only 
recently at the microscopic level in disordered systems 
[ill Hq l . Second, the shear stress perturbation alluded 
to in (iv) is computed exactly within the framework of 
tensorial linear elasticity for an isotropic incompressible 
material as reported in Ref. This provides the ex- 

plicit Green function, Gxzxz, relating the stress variation, 
5a, at any point in the system, to the xz component of 
the plastic strain e^\{x' , z'}; t), associated with the plas- 
tic event localized at {x', z'}. Using the simpler notation 
G for this function yields 

5a{{x,z];t) = 2^i j dr' G{x,x' , z, z'y\{x' , z'];t) (1) 

The shear modulus /i has been exhibited for conve- 
nience. In a 2D infinite system, G decreases as G(r) = 
l/7rr^ cos(40) - in cylindrical coordinates {r, - (in 
agreement with refs [ly, Lj{ ) . In general, its precise form 
depends on the specific geometry of the system : infinite, 
periodic or confined between two rigid walls Sum- 
ming up at this point, the evolution of the shear stress 
field results from the global elastic loading 7 plus the 
perturbations induced by the localized plastic events: 

dt(T{{x,z],t)= pL^ + 2^i j dr' G{x,x\z,z'y\{x\z'};t) 

(2) 

The last part of the modelization is the choice of a dy- 
namical law for the plastic events, i. e. the feed-back law 
relating the plastic relaxation eP^{{x, z}]t) to the stress 
field a{{x\ z'}\t' < t). As in the Princen model, we 
choose a local relation with a threshold stress value ay ■ 
In addition, an intrinsic time scale r is introduced to 
describe the dynamics of the event. We anticipate that 
this will lead to a shear rate dependence of the dynam- 
ical structure in the flow, driving the system away from 
the critical quasi-static limit (self-organized criticality in 
a related quasi-static model was reported in to a 

more homogenous situation at large shear rate. Another 
important outcome is that the local stress may exceed the 
yield stress ay for a finite time interval so that the aver- 
aged stress can also grow beyond this value, as observed 
experimentally. 

There are actually many possibilities to introduce such 
an intrinsic time scale for plastic events, and few guides 
as to how should do so. We make here a simple arbi- 
trary choice and assume that the system locally alter- 
nates between a purely elastic state and a plastic state 
(during which stress is released), with finite transition 
rates : t^^^^j is the rate of transitions from elastic to 
plastic, while the reverse transition is characterized by 
a time Teiast- Since plastic events only occur above the 
yield stress ay, we take Tpiasticr) = 00 if locally a<ay. 
We otherwise assume for sake of simplicity fixed values 



for the Teiast and Tpiast, independent of the local stress. In 
order to finalize our model, we eventually have to quan- 
tify the amount of plastic strain released in an event and 
simply assume a Maxwell, visco-elastic like relaxation of 
the material in the plastic state ipiast — 2^''^' with r a 
mechanical relaxation time. All the previous discussion 
is best summarized by introducting a "state variable" 
n{x,z) such that n = 0/1 identifies the elastic/plastic 
state : 

eP^{{x,z},t) ^ ^n{{x,z},t) a{{x,z},t) 
n{{x,z},t) : 1 O'^'l (3) 

if (T>(Ty Vcr 

Equations (O and ((SJ constitutes our minimal starting 
point to describe the dynamics of yield stress materi- 
als under flow. Note that, as in the somewhat related 
analysis of Langer neither the stress nor the state 
variable are convected by the displacement field within 
the present simplified model. 

Before turning to their resolution, eqs (|2I3|I are made 
dimensionless using ay and r as stress and time units. 
An important point emerging from this procedure is that 
the shear rate only appears in the form of the ratio j /jc, 
with 7c = ay/fiT. In this dimensionless form, our model 
therefore therefore points out to a very general scenario, 
in which specific microscopic details are embeded in the 
precise values of cry and 7c, as already suggested by some 
experiments |5Ll2ll|. 

The dynamical equations in this dimensionless form 
have been solved numerically, by discretizing the mate- 
rial into blocks of elementary size a. A pseudospectral 
method is used, which allows us to express easily the 
stress increments in reciprocal space at each timestep. 
On the other hand, the state variable n(i a,j a) in the 
block {i, j} evolves in real space according to the stochas- 
tic laws enonced above. We have focused on two ge- 
ometries of iV = (L/a)^ blocks : a biperiodic geom- 
etry and a confined one where the system is bounded 
by two rigid parallel walls. Practically we have chosen 
Tpiast = Teiast — T for the results reported here. 

We first quote the results for the biperiodic system. 
In Fig. n we plot the results for the macroscopic flow 
curve, which displays the essential features observed in 
experiments. First a plateau is found at small shear rate, 
defining a macroscopic yield stress at vanishing shear 
rates. The latter is found to be lower than the micro- 
scopic yield stress ay, and also lower than the related 
peak value of the stress versus time at small shear rates 
(see inset in Fig. . A different regime is found at large 
shear rates, where a Newtonian behaviour is recovered. 
The dynamics of the time dependent shear stress is also 
quite different in these two regimes. In particular, rel- 
ative stress fluctuations around the mean value increase 
as the shear rate decreases (not shown), in agreement 
with observations in experiments and simulations 0, 0] ■ 
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FIG. 1: Shear stress versus shear rate in units of ay and 
7c = (Jy I \i.T (log- log plot). The various symbols correspond to 
four different system sizes (from 4 x 4 to 32 x 32 blocks) . Inset; 
time dependent stress at the low shear rate corresponding to 
the arrow. The dashed line corresponds to the microscopic 
yield stress ay ■ 



A zoom on the dynamics at shorter time scale actually 
shows that at small shear rate the stress exhibits succes- 
sive periods of elastic raise and abrupt drops, as observed 
in the quasistatic limit in various systems 0|. These 
drops encompass many events, constituting "bursts" of 
correlated plastic activity. More interestingly these dy- 
namically correlated events are also highly correlated in 
space, as emphasized in Fig. |21where the spatial distribu- 
tion of the cumulated plastic activity is plotted for a given 
succession of plastic events. This figure clearly shows 
that while at high shear rate plastic events are spatially 
decorrelated, a correlation pattern shows up as the shear 
rate is decreased, leading to the development of long-lived 
"fragile" zones in the system where nearly all the plastic 
activity takes place. This graph therefore suggests the 




■1 JU 



FIG. 2: Spatial distribution in {x,z} plane of the cumulated 
plastic activity for a A'' = 16 x 16 system after 21^ = 512 
plastic events. From left to right the dimensionless shear rate, 
7/7C, is : 5.10"^*, 5.10~^, 0.05, 2. Gray levels correspond to the 
number of plastic events. Arrows indicate the imposed shear 
direction. 



development of a shear rate dependent length scale in 
the system, which grows at small shear rates. In order 
to get more insight into this aspect, a possible route is to 
measure the length via measurements of correlation func- 



tion. This is however a difficult task in general [22| and 
we have followed a different strategy here, analogous to 
finite size scaling. Namely, since such a length is associ- 
ated with the correlation of plastic events during a macro- 
scopic stress drop, it should show up in the statistics of 
stress drops. To this end, we have computed the mean 
stress drop in plastic events. Act, as a function of shear 
rate [i^l- Results are shown in Fig. 01 for various system 
sizes. Let us first discuss the inset which exhibits the bare 
results for the average stress drop normalized by the av- 
erage stress. Act = Act(7)/ct(7). Three different regimes 
can be identified for all system sizes : two plateaus at 
large and small shear rates and an intermediate regime 
relating these two. We remark that the transition be- 
tween the intermediate and "saturation" regime at low 
7 shifts to lower shear rates when the size of the system 
is increased. The system size dependence of this tran- 
sition at low shear rate is best evidenced using a rescal- 
ing procedure, in which all different curves are found to 
coincide in the variables {A^7/7c, Act(7)/Act(0)}. This 
rescaling in the 7 variable actually allows us to quantify 
the correlation spatial effect observed above. Assuming 
the existence of a shear rate dependent correlation length 
^(7) in the system, a saturation effect is expected for 
the mean stress drop when ^(7) reaches the system size, 
iV^/^a. The rescaled graph indicates that such a satu- 
ration occurs for a fixed value of N'j/'^c which suggests 
^(7) ~ 7~", with a ~ 1/2 from these data. Our model 
therefore explicitly yields indication of (at least) one di- 
verging length scale at small shear rates, a feature to our 
knowledge absent in previous studies of yield stress fluids. 
Interestingly, the transition between the intermediate to 
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FIG. 3: Mean stress drop (normalized by the average stress), 
Act, as a function of the renormalized dimensionless shear 
rate, N^f-yc- Results in the inset are shown for five differ- 
ent sizes (from 2 x 2 to 32 x 32 blocks from top to bottom), 
while in the main graph all curves are rescaled using vari- 
ables {Nj/'jc, A(t(7)/A(t(0)} to emphasize the scaling of the 
transition between the small shear rate plateau and the in- 
termediate regime. The dotted lines in the inset sketch the 
separation between the three dynamical regimes. 

the large shear rate regime on Fig. 13 (right dotted line in 
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the inset) occurs roughly at the characteristic shear rate 
7c, independent of system size, as for the macroscopic 
flow curve in Fig. ^ From these first results, our model 
yields a flow behaviour with three different regimes as 
sketched on Fig. ^ (i) for 7 > 7c (or cr > try), the 
blocks uncorrelated in their dynamics and the flow is ho- 
mogeneous ; (ii) for 7 < 7c, correlations increase up to 
a correlation length ^(7) which diverges algebraically at 
small shear rates; (iii) at very low shear rates, the corre- 
lation length saturates at the size of the system, leading 
to a quasi-static dynamical behaviour. 
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where the deformation of the system takes place. These 
bursts take place preferentially close to the walls. At 
high shear rates, this correlation length is comparable to 
the size of the individual elements which flow indepen- 
dently from one another. These features are essentially 
compatible with recent observations in experimental or 
numerical systems : localization of the time-averaged de- 
formation 13,0,0, ^ IE IE 1^ 1 intermittency at low shear 
rate 0,IEl3,0|) a diver ging length scale at small shear 
rate in granular systems 24]. Moreover numerical simu- 
lation of glassy systems j^J show that flow heterogeneities 
occur for global shear rates such that a < ay, a conclu- 
sion which is recovered within our minimal model. Al- 
though our model should be refined to take into account 
convection and the full tensorial nature of the problem, 
the present early results suggest that the generic behav- 
iors observed in the experiments and molecular simula- 
tions originate in a minimal number of ingredients. This 
opens the possibility for a coherent and robust scenario 
for the slow flow behaviour of disordered materials. 



FIG. 4: Sketch of the emerging flow scenario in the 
{'y/'jcyi/N) plane. Successive transitions from a homoge- 
neous flow to an organized and a finite size regime occur as the 
correlation length ^ grows from the block size to the system 
size, as the shear rate decreases. 

We have also studied a confined geometry where two 
rigid walls bound the system in the z direction. A deli- 
cate technical point is then the calculation of the Green 
function, which shows that shear stress perturbation is 
amplified close to the walls ^E|- Essentially, the picture 
in the confined geometry is very similar to that of the 
biperiodic system described above (Figs 1, 3, 4). One 
important specific feature however concerns the localiza- 
tion of the flow : while at high shear rate the flow is 
homogeneous, at low shear rates the plastic bursts occur 
preferentially close to the walls. In this last regime, the 
average flow corresponds to an increased shear rate close 
to the walls, but this "localization on average" of the 
flow is only part of a complex spatio-temporal pattern. 
A more detailed analysis of this regime is left for a future 
publication. 

To sum up, we have proposed an athermal elasto- 
plastic model for the flow of yield stress systems, con- 
structed on the basis of two generic ingredients: localized 
plastic events, occurring above a microscopic yield stress 
with a finite duration, and an otherwise elastic behaviour 
of the material (including redistribution of stress during 
the events). These two ingredients lead to a complex 
spatio-temporal behaviour of the system at small shear 
rates. More precisely a correlation length is exhibited 
which diverges at small shear rates, corresponding to in- 
termittent collective events (correlated bursts of plastic 
events) , leading to the creation of (long live) fragile zones 
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